arXiv: 1505.02470vl [quant-ph] 11 May 2015 


i 

Coherent Quantum Control of S2 -H- Si Internal Conversion in Pyrazine 
via So —> S2/S1 Weak Field Excitation 

Timur Grinev [1] 

Department Of Chemistry, Chemical Physics Theory Group, and Center for Quantum Information and Quantum Control, 
University of Toronto, Toronto, Ontario M5S 3H6, Canada 

Moshe Shapiro 

Department of Chemistry and Department of Physics, University of British Columbia, Vancouver, British Columbia 
V6T 1Z1, Canada, and Department of Chemical Physics, Weizmann Institute of Science, Rehovot 76100, Israel 

Paul Brumer 

Department Of Chemistry, Chemical Physics Theory Group, and Center for Quantum Information and Quantum Control, 
University of Toronto, Toronto, Ontario M5S 3H6, Canada 

Abstract 

Coherent control of internal conversion (IC) between the first (Si) and second (S 2 ) singlet excited electronic 
states in pyrazine, where the S 2 state is populated from the ground singlet electronic state Sq by weak field 
excitation, is examined. Control is implemented by shaping the laser which excites S 2 . Excitation and 
IC are considered simultaneously, using the recently introduced resonance-based control approach. Highly 
successful control is achieved by optimizing both the amplitude and phase profiles of the laser spectrum. 

The dependence of control on the properties of resonances in S 2 is demonstrated. 


I. INTRODUCTION 


Coherent quantum control 2|, 3] has been extensively studied for a wide variety of systems and proven to be 


a useful approach to controlling properties of atomic and molecular systems. For example, in bound systems it has 
been used to suppress spontaneous emission from a manifold of states [d], and to control radiationless transitions in 
collinear carbonyl sulfide OCS [fj| and in pyrazine C4H4N2 

fifl 


Christopher et al. examined |]6j,|8| radiationless transitions in pyrazine from the S 2 to the Si electronic state and 
controlled the process by optimizing the superposition states belonging to S^. The problem was first studied [6] using 
a simplified four-mode model for the pyrazine vibrational motion |9|. The optimization technique used showed the 
possibility of performing active phase control of S 2 O Si interconversion, and that this control is directly related to the 


presence of overlapping resonances 


motion of pyrazine 


[lO, 11] i 




in the S 2 manifold. Subsequently [7, 8], the full 24-dinrensional vibrational 


12] was considered, and the dynamical problem solved using an efficient Lowdin-Feshbacli QP- 
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partitioning approach. Previous control results were fully confirmed and refined, proving the high controllability of 
S 2 -O- Si internal conversion by actively exploiting the effect of quantum interferences which was shown to rely on the 
presence of overlapping resonances. 


In Refs. ( 6 , | 8 | coherent control was implemented for pyrazine that was already prepared in the excited S 2 
state. The Sq -4 S 2 excitation process was not considered, assuming instead that the excited states in S 2 were already 
populated. Recently, we showed the possibility of performing effective coherent control in a simple IBr diatomic model, 
where we explicitly included the exciting laser in an approach that simultaneously considered excitation and decay 
to a continuum [13L 14]. In that case we introduced an optimization schemes different from the simple one used in 


Refs. |] 6 , 81 and demonstrated the reliance of control on overlapping resonances. Below we considerably generalize 
this study to pyrazine, explicitly introducing the laser to excite the 24-dimensional Si + S 2 vibronic pyrazine model 




and using the same control and optimization schemes as for IBr 
controllability on the properties of the S 2 resonances in pyrazine. 


43], Signi 


Significantly, we confirm the dependence of 


This paper is organized as follows. Section |TT] reviews the theory explicitly accounting for the exciting laser 
in the weak field limit. Section m introduces the coherent control approach for the S 2 population, points out its 
connection with the properties of S 2 resonances, and provides additional details of the approach. Section HVl provides 
computational results for control of pyrazine internal conversion. Section [V] provides a summary and conclusions. 


II. So -4 S 2 EXCITATION AND S 2 <4 Si INTERNAL CONVERSION 


Below, |«) denotes vibrational states belonging to the S 2 electronic state, with corresponding projection operator 
Q = |«)(«|- Since the |k) states are not eigenstates of the full Pyrazine Hamiltonian, the system evolves in time 

if it were prepared in these states. Hence, such states are termed resonances. The states |/3) denote vibrational states 
belonging to the Si electronic state, with P = |/3) (/3| being the associated projection operator. The full vibronic 

states, which are eigenstates of the full Pyrazine system, are denoted I 7 ), so that P + Q = I = ^ |t)(tI- 


A. Time Evolution of the System Assumed Already Excited 


In Refs. 


Sq -4 S 2 laser excitation is assumed to allow preparation of a superposition of |re) resonances: 


l*(o )> = X>Ik'>- 

k' 


( 1 ) 


The dynamics of internal conversion was then described by an action of the propagator U(t) on |T(0)): T(t)) = 

t/(t)|4'(0)). Because |y) are exact states of the system Hamiltonian, the spectral resolution of the evolution operator 












3 


U(t) = exp(—iHt/h) is U(t) = exp(— iE 1 t/h)\'y){y\. This gives 

I V(t)) = E c "'E exp ( - *^V ? 0<7|K , }l7} ='^2a 1 exp{-iE 1 t/h)\'y), (2) 

k,' 7 7 

where a 7 = X) K , c K /( 7 |/c'). 

The S 2 electronic state population Pg 2 at time t is an observable defined by the projection operator Q onto the 
state |T(t)}: 

p s 2 it) = (y(t)\Q\V(t)) = E a*,(t)ay>(t)Qy l7 », (3) 

7 , j7 // 

where a 7 (t) = a 7 exp(— iE^t/K) and Qy t y' = (7 / |Q|7 // )- Equation Q can be rewritten in matrix form as: 

Ps 2 (t) = ati^Qi^a, (4) 

where a is a vector with a 7 components, e ±lEt / n are square diagonal matrices composed of exp(±iE-yt/h) values, and 
Q is a square matrix with Qy y< matrix elements. 

Since Q = l K )( K l> the matrix elements Qy n » = {l'\QW) = ]Ck(7 , | k )( 7 ,, | k )*■ Introducing the matrix R 
with R 7iK = ( 7 |k), then Q 7 ' >7 " = J2 K Ry^R*",k = = [RR%>,y, giving 

Q = RR t . (5) 

In turn, according to Eq. ©, the vector a can be written as 

a = Rc, (6) 

where c is a vector composed of c K > coefficients. Inserting Eqs. © and (f6j) into Eq. © gives: 

Ps 2 {t) = clRt^^/^RRtg-^/^Rc = c t M c t(t)M c (t)e = c^K c (t)c 

= E c >*"K c K ’, K ''(t) = E M 2 K° K , tK ,{t) + E <Z'C*»K 0 K ' tK „(t), (7) 

K,' ,K," K,' k'^k" 

where M c (f) and K c (f) matrices are defined as 

M c (f) = RV iBt / R R, K c (t) e M c t(f)M c (t) = rV^/'RRV^^R. (8) 

The matrix elements of M c (t) have the form 

= E<^l7>(7l« / ) exp (-*^i/^) = («l U(t)W), (9) 

7 

being matrix elements of the U(t) propagator operating between the resonances |k) and \k'). According to Eq. ©, 
M £ K ,(t ) ^ 0 for k ^ k', only if there is at least one state |y) such that (k| 7 ) ^ 0 and ( 7 |k') ^ 0. If so, then resonances 
\k) and \k') are said to be overlapping. This resonance overlap property is crucial for M c (f) nondiagonality which, in 
turn, provides K c (f) nondiagonality, which allows efficient phase control of Ps 2 (t ) in Eq. © by means of phases <p K i 
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of complex coefficients c K / = \c K >\ exp(i(p K ') 


Such phase control is termed active control, in contrast to passive 


control, which is control via the |c K '| amplitudes only. 


In the case of pyrazine, which has 24 vibrational degrees of freedom, there is a large number of |y) states 


M 


To make the computations feasible, instead of exact states, a set of approximate coarse-grained states is used 
to compute the time evolution. Specifically, the energy axis is divided into small bins I a , of size A a , center energy E a 
and density of states p a . The projector onto the coarse-grained state |a) is then defined as: 


|q:)(q!| — (l/(Pa^a)) ^ ' |7)<7|, hence \Jp a |ck)( ck| yjp a A a = 51 ItMI- 

7 e/c 7 e/ a 

Thus, the coarse-grained state |a) effectively replaces all the |y) states in the bin I a . Numerically, the weighted 
states |a) = y/ p a A a \a) and their overlaps with resonances \k) are available through our iterative solution method for 
pyrazine, based on QP-partitioning algorithm (described in detail in Ref. j^]), giving 


|a)H = 55 ItXtI- 

7 £l a 

All the 1 7 ) states belonging to the same bin I a are treated as one effective state |a); so that 


( 10 ) 


= 55<k|7X7|*'> exvi-iE-tt/K) = 55 55 ( K l7)(7K) expi-iE^t/h) 

7 cx. 7 G/q 

Ki'^2(n\a)(a\K , )p a A a ■ — exp {-iE y t/K). 

a A “ 7 ei a Pa 


( 11 ) 


The remaining inner sum over 7 £ I„ in Eq. CD is approximated by a corresponding integral: 


7 7 7 rE a + A q /2 

—— —exp(— iEjt/h) « —— / dEry exp(— iE^t/h) 

_ , Pa , 

7 


/ E a — A a /2 


, . _ s sm(A a t/(2h)) , . 
= exp(-,E„t/B) EUX - a r„(<), 


( 12 ) 


giving the final coarse-grained expression for M f (: K , (f): 


« 55( K l a X a l K ')T-«(X = (k| 

a 


55 r«(t)|a>(a| 


!«')■ 


(13) 


The quantity in the square brackets is the coarse-grained approximation to the 17(f) propagator, and the sum is 
over all available |cf) states. Equation m is accurate for the evolution times which are not too large, i.e., when 
\ T a(t)\ = \sm(A a t/(2h))/(A a t/(2h))\ ~ 1 , implying that \t\ -C 2h/A a . The resonance overlap phenomenon and the 
need for nonzero coarse-grained off-diagonal M K ,(t) discussed above remains the same, except that the |y) states are 
replaced by |a) states. 
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B. Time Evolution Due to Laser Excitation 


Consider now the result of single photon excitation from the ground electronic state So, which produces the 
excited time-dependent wavepacket, as a superposition of I 7 ) states (here the subscript p denotes pulse): 

I’M*)) = ^2,b 1 (t)eyLp(-iE 1 t/K)\i), (14) 

7 

where & 7 (f) coefficients are, in general, time-dependent. 

The S 2 electronic state population at time t is given by: 


PsAt) = <M*)IQIM*)> = ~ h i '(*) 

7 ', 7 " 

where & 7 (f) = 6 7 (f) e~xjp(—iE~/t/h). Equation (11511 can be written in matrix form as 


(15) 


Ps 2 (t) = h\t)g iEt / h Qg- lEt/n b(t), 


(16) 


where b(t) is a vector composed of & 7 (f) components. 

If the exciting laser pulse is weak, first-order time-dependent perturbation theory is applicable, and the 6 7 (f) 
expansion coefficients in Eq. m can be written as 


6 7 (£) = (i/K)('y\n\g)e p (U'y,g,t), (17) 

where p is the dipole operator, | g) is the ground vibrational state on So, w 7iS = ( E 1 — E g )/U, and £ p (oj ltg ,t) is the 
finite-time Fourier transform of the e p {t ): 

e p (u> 1 } g,t) = / dt r £ p (t') exp(^a; 7J^ ^ , ). (18) 

J — OO 

Eq. (1171) can be written in matrix-vector form as 


b(t)=^£{t), (19) 

where p is a square diagonal matrix composed of (i/h)(”/\/j,\g) values, and e(£) is a vector composed of £ p {uj ltg ,t) 
components. 

Inserting Eqs. ([5|) and (1191) into Eq. © gives, for the Ps 2 ( t ) population, 


PsAt ) = A (t) A §. iEt ^ n PPA e~ iEt / n n e{t) = A{t)M E At)M%t)£(t) = A (t)K s (%(£) 


^ ) l^p( t * , 7',g’ *) I *^7',7'(*)^ _ e p( w 7',3’*) £ p( W 7",S’*)*^7',7" (*)’ 

7' 7'^7" 

where M e (f) and K E (£) matrices are defined as 

M E (t) = Rt e~ iEt l n pL, K E (£) = M Et it) M E (t). 


( 20 ) 


( 21 ) 
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Since p and e ±lEt / h are diagonal, the only source of nondiagonality in Eqs. and (EB for K £ (t) is Q = RRt. 
Thus, phase control via the phases of complex e p (u; 7ig , t) = |e p (w 7>g , t)\ exp(i</> 7 (t)), depends solely on properties 

of Q. 

A few comments are in order. First, R is a rectangular matrix, with each K th column composed of overlaps 
i? 7)K = ( 7 |ft) of the resonance |k) with all |y) states. On the one hand, each resonance, being broadened in energy, 
has more than one nonzero ( 7 |k) term in its Kth own column. On the other hand, if resonances |k) and |k/) overlap, 
then they have at least one common |y) such that, for this |y), both i? 7iK 0 and i? 7 , K ' ^ 0 simultaneously. 

Second, all nonzero ( 7 |ft) components of each column in the R matrix that are related to one particular 
resonance |k) form a square block centered along the main diagonal in the resulting Q = RRf matrix, filled by terms 
<3 7 ' )7 " = ( 7 '|av}(k| 7 "). Thus, Q displays block-diagonal structure. Since each block dimensionality is larger than one 
due to resonance energy broadening, nondiagonal matrix elements in these blocks are generally nonzero, contributing 
to K e (t) nondiagonality, and thereby providing Ps 2 (t) phase control associated with the energy broadening of each 
particular resonance. This kind of control will be discussed below. Furthermore, if resonances |k) and \k') overlap, 
then the corresponding blocks overlap, so that the Q matrix acquires a non-block-diagonal structure. In this case 
Q 7 / i7 // matrix elements belonging to two blocks simultaneously are a sum of terms borrowed from each block (produced 
by its corresponding resonance): Q 1 \ 1 " = ( 7 , |k)(k| 7 ,/ ) + Similarly, in the case of overlap of N blocks, 

the sum contains N terms: (7 , | K }( K 'l7 ,, )• As will be discussed below, the resonance overlap effect 

greatly increases the overall phase controllability in comparison with a pure resonance energy broadening effect. 

The nondiagonality in this section (see above), is very different from that discussed in Sect. Ill Al Specifi¬ 
cally, in Eq. 0 - for the case when the system is already assumed to be excited, control is performed by means 
of the c K > coeeficients, so that a = Rc, giving K c (t) = Rfe I ' Bt / ri RRle _zEt / n R [Eq. (0)]. This greatly simplifies 
the K c (f) nondiagonality dependence, effectively removing the resonance broadening effect and leaving only reso¬ 
nance overlap as the crucial effect that provides nondiagonality, i.e., phase control. By contrast, in this section, 
K £ (t) = p^e zEt ^ h IllOe~ lEt ^ n p [Eq. (T2T1) ] and nondiagonality is provided only by the Q = RRi matrix itself, whose 
nondiagonality, responsible for phase control, depends on both resonance broadening and resonance overlap effects. 


It can be noted that M £ (t)e(t) in Eq. (l20l) is a vector composed of components 


(«|7> exp(-*E 7 t/^)-(7|/r| 5 ) 


( 22 ) 


(K'l’^pW) — ^ ^ S p (uj 7;3 , f)Af K . 7 (t) — ^ 7,gi 

7 7 

In the case of pyrazine, transition dipole matrix elements for the Sq —> <Si excitation are an order of magnitude smaller 
than for the So —> S 2 excitation Q- 8 , 15 1 , thus allowing the following “doorway” approximation: 


< 7 |mI 9) = (7 1 (P + Q)t<\g) = g) ~ ^(i^^hAg)- 

f3 K, K, 


(23) 


Equation (1231) indicates that the excitation to a full vibronic state [ 7 ) takes place by means of a preliminary intermediate 
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transition to a manifold of |re) resonances. Inserting Eq. (l23l) into Eq. (l22l) gives 


7 

which can be rewritten as 


(«|7> exp (-iE^t/h)^ 




53 £ p( w 7,s^)( K l7){7|K , ) ex P(-*- B 77 ?i ) 


(24) 


(25) 


In order to make the computations below feasible, we introduce here a coarse-graining procedure for the quantity 
in square brackets in Eq. (H>D- This procedure is similar to the one made in Ref. {tJ, taking into account Eqs. m 
and m- Namely, is written as J2 a X) 7 e : 


F y, eyiv(-iE^t/%) tt'yE p {u <Xt g,t){K\a){a\K')p a A. ol ■ — y — exp(-*E 7 t/?i) 

a 7e/o, a A “ 7 e/ Q 

~ X] £ p( w “. 9’ t)(K|a)(a|K')Ta(*)? 

Q! 

where = (E a — E g )/Ti. Inserting Eq. (1^1) into Eq. K5l) gives: 


(«:|4'p(t)) ss 'y ' E p (uj at g, t) 


K\a) r a (t) ^ y{a\K'){n'\n\g) 


— Ep(u}g,g, t)M^ a (t). 


Below, a superscript a indicates the coarse-grained nature of the corresponding values. Here, the quantity 


y tW\v\9)W) 


(26) 


(27) 


(28) 


= («|a) T a(t) - y(a\K)(K\p\g) = (k\ [r„(t) \a){a\} 

K' 

is a coarse-grained version of [Eq. (p?2]) ] and depends only on the material system properties. If one defines 

\ y(aW)W\9\g) = ^\g\g), K, k = (a|«), (29) 

k' 

then 


M e ’“(t) = R at ^ a (t)/i a , K e ’ a (f) =M e ' at (t)M e ’ a (f), (30) 

where r a (t) is a square diagonal matrix composed of r a {t) values, and is a square diagonal matrix composed of 
(i/h)(a\g\g) values. Then the Ps 2 {t) population in terms of coarse-grained values becomes 

Ps 2 (t) =£ a \t)M e ’ a \t)M e ’ a {t)£ a {t) = E a \t)K e ' a (t)E a {t) 

= y\E p (u} a ',g,t)\ 2 K e a ? a ,(t)+ y El(u a , jg ,t)E p (u a ,', g ,t)K%“ a ,,{t), (31) 

a' 

where £ a (t) is a vector composed of s p (u> a .g, t) components. 

The quantities g a and r a (t) are diagonal matrices, so the only origin of nondiagonality in Eq. (1301) for K e, “(t) 
and Eq. ([3D is via the Q“ = R Q R°^ matrix, composed of Qf, a „ = {a'\Q\a") matrix elements. Hence, all the Ps 2 (t ) 
phase control considerations from above remain the same, except that |y) states are replaced by |a) states. Namely, 
phase control is driven both by resonance energy broadening and resonance overlap. The resonance overlap effect, 
providing a non-block-diagional structure of Q“ and K e,a (t), strongly enhances the effect of resonance broadening. 










III. COHERENT CONTROL OF PYRAZINE INTERNAL CONVERSION 


Section III Bl above describes resonance broadening and resonance overlap, two effects related to Q (Q“) and 
K E (f) (K e, “(f)) nondiagonality. Here, a control scheme based on resonance broadening is discussed in Sect. IIII Al 
Section UlI Bl discusses a control scheme relying on presence of resonance overlap. 


A. Control Associated with Single Resonance 


In the case of pure resonance broadening without resonance overlap, one particular resonance |/c) has nonzero 
R-y,n = ( 7 |k) terms for some specific set { 7 }^ of |y) states. This results in the simplified expressions for K e (f) matrix 
elements for this { 7 } K set, with the summation over k reduced to a single term 


= |M* 7 ,(t)| 2 , A-, > 7 "(t) = M£ 7 ,(i)M* 7 „(t) 

for the diagonal and nondiagonal matrix elements, respectively. 


(32) 


The probability Ps 2 {t) [Eq. (EOl) ] is a quadratic form of complex time-dependent variables £p(w 7 i 9 ,f). When 
the pulse is already over (at t = T over ), these values become infinite-time Fourier transforms of this laser pulse at 
different frequencies, e p (u ; 7iS ); they are no longer time-dependent for t > T over . Here we use the so-called absolute 
control scheme for Ps 2 (t) optimization, with the K e (f) matrix given in Eq. (I.32f) . Namely, Ps 2 (t ) is optimized at a 
desired optimization time t = T, while keeping the total energy of the pulse at 2ttEq: 

|£p(w 7j g)| 2 = £ f £ = 2itE 0 . (33) 

{7}* 

This is done by introducing the corresponding Lagrange multiplier A " 4 (superscript A denotes absolute) with the 


corresponding optimization function at time T defined as: 


Ps 2 A {T, e) = e t K E (T)e - A A (£ f £ - 2nE 0 ). 


(34) 


We then search for /Av ' 4 (T, e) extrema with respect to e: 


dPs 2 A (T,s) 
d Re [£p(w 7 , s )] 
dPsf(T,e) 


= 0 , 


= 0, 7 = 1,...,TV {7}k 


(35) 


<9Im[£ p (w 7>g )] 

where V{ 7 i is the number of I 7 ) states in the set { 7 }^. Conditions in Eq. (1551) . applied to Eq. (1M1) . lead directly to 
an eigenvalue problem 


K E (T)e = \ A £. 


(36) 


which provides a set of eigenvalues A ' 4 and corresponding eigenvectors e with a unit norm (e'I'e = 1). Multiplication 
of these e eigenvectors by \J2-kEq provides the required optimized solutions. 
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The K e (t) matrix is such that all but one of its 1 Vj 7 \ eigenvalues are equal exactly to 0, while its last eigenvalue 
is equal to the sum of its diagonal elements: 

A^ = 0, n = 1,... — 1; = £ \M^(t)? ■ (37) 

{7}* {7}- 

This is an analytical property of the K £ (t) matrix in Eq. (13211 . so that a numerical solution of the eigenproblem in 
[Eq. ((5511 ] is not required. Specifically, for any time T, Ps 2 (T) can be set to zero, using the eigenvector corresponding 
to zero eigenvalue. In terms of the coarse-grained |a) states, the results are the same with the { 7 }^ set replaced by 
{«}«■ 

Given the simplistic nature of this solution, numerical results are neither necessary nor are they provided below. 
Note, however, that this type of control is possible only if the system displays isolated resonances. This can be 
the case in small molecules; large molecules such as pyrazine, however, display overlapping resonances throughout 
the spectrum, with highly unlikely regions of isolated resonance. Such systems can be controlled via an alternate 
mechanism, discussed below. 


B. Control Associated with Overlapping Resonances 


Here we consider a second different control scheme, termed relative control. Namely, we optimize the ratio of 
Ps 2 (t) populations at times T 2 and Ti, where T 2 > Tf > T over : 


,fi P S 2 ( T 2 ) , 

A =-—— —> max, mm 


(38) 


PsAT,) 

(where superscript R denotes relative). One can optimize the value of Ps 2 (T 2 ), keeping the value of Ps 2 (Ti) constant 
5 and equal to some predefined value Pq. Here fixed Ps 2 (Pi) = Po assures that enhanced (or diminished) Ps 2 at the 
target final time T 2 does not simply result from a stronger (or weaker) field that simply achieves control by affecting 
the amount of S 2 excited. To do so, we consider the optimization function 


P* r (T 2 ,T i,£) = Akat 2 )£-x r (Ak s (t 1 )£~Po), 


(39) 


where \ R is a yet unknown Lagrange multiplier. We then find P s ' 2 (T 2 , Ti, e) extrema with respect to e leading directly 
to a generalized eigenvalue problem: 


K e (T 2 )£ = A fl K £ (Ti)e. 

Multiplying Eq. (HOI) by £' from the left gives 


(40) 


e t K £ (T 2 )e = AVK e (Ti)e. 


(41) 
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The X R is real and positive because £ 1 'K e (T 2 )£ = Ps 2 (T 2 ) and £ t K E (T'i)£ = Ps 2 (7i) are real positive values. Dividing 
Eq. (HTl) by £^K e (Ti)£, yields X R = Ps 2 {T 2 )/ Ps 2 {Ti), i.e., X R is the optimized ratio of the populations of interest [Eq. 

mi 

The K E (t) matrix determinant is generally nonzero at every time t, so that K E (t) always has an inverse [K £ (f)] _1 . 
This allows transformation of the generalized eigenvalue problem in Eq. ( 1401 ) into an ordinary eigenvalue problem. To 
do this, we multiply the left and right sides of Eq. fiUl) by [K e (Ti)] 1 from the left: 

R e (T 2 , Xi )e = X R e, (42) 

= [K e (Ti)] _1 K e (T 2 ). (43) 


The solution to the eigenproblem in Eq. (H^l) for times T 2 > T) > T over is dependent only on the properties of the 
material system. Moreover, this solution is the best possible in the weak field case, i.e., it is optimal 17]. Specifically, 


the maximal and minimal eigenvalues X R provide the entire achievable range of Ps 2 (T 2 ) / Ps 2 {T{] for a given T 2 and 
Ti, obtained using the corresponding eigenvectors e. 


In terms of coarse-grained states |a), e is replaced by £ a , and K e (t) is replaced by K e ’“(t), giving the following 
coarse-grained version of the optimization problem: 

R E ’ a (T 2 ,Ti)£ a = X R ’ a £ a , (44) 

R e ’“(T 2 ,Ti) = [K £,a (Ti)] _1 K E ’“(T 2 ). (45) 


In addressing this problem computationally, we encountered numerical instability in Eq. (1441) if the number of 
|a) states is relatively large (150-180). Namely, the condition number of K E,a (t) tends to become very large, resulting 
in an ill-conditioned matrix, preventing accurate numerical construction of R E,a (T 2 ,T 1 ) [Eq. (1451) ] and its subsequent 
diagonalization. To overcome this problem, we partitioned the energy axis into a limited number of Na bins in Eq. 
m, as discussed in the Appendix, giving further broadened |A) states. 


Using these further broadened |A) states allows us to reformulate the eigenproblem in Eq. (l44l) as 


R e,A (T 2 , Ti)£ A = X R ' A e A , 

(46) 

R E ’ A (P 2 , T!) = ]K E ’ A (T!)] - 1 K E ’ A (P 2 ), 

(47) 


where the states |a) in Eqs. m and (l45l) are replaced by the further broadened states | A), as described in the 
Appendix. 


C. Numerical Correlation between Controllability and Resonance Overlap 

In general, effects of resonance energy broadening and resonance overlap are mixed together in the structure 
of the Q a and K E,A (f) matrices. To quantitatively estimate the K E,A (t) nondiagonality, providing phase control, we 
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utilize the Hadamard measure: 


H (K E ’ A (f)) = det (K E ’ A (f)) /det (diag (K E ’ A (f))), (48) 

where det denotes a determinant, and diag is the diagonal part of a matrix. Thus, det (diag (K E,A (f))) = 
Since K E,A (f) is a Hermitian positive-definite matrix, both det (K E,A (f)) and det (diag (K E,A (f))) 
are real and positive. Furthermore, det (K E,A (f)) < det (diag (K E,A (f))), giving 

0 < H(K £ ' A (t)) < 1, (49) 


where the equality applies if and only if K E,A (f) is strictly diagonal. 
The determinant of R e ’ A (T 2 ,Ti) can be expressed as 


det (R e > a (T 2 ,Ti)) = det [K e ’ a (T 0] 1 K E ’ A (T 2 ) = det (K e ’ A (T 2 )) /det (K e - a (Ti)) . 


(50) 


Hadamard-likc measures of non-diagonality for R E,A (f) are introduced in a similar manner: 


Hr (R e,A (T 2 ,Ti)) 
He (R E,A (T 2 ,Ti)) 


_ det (R E ’ A (T 2 , Ti)) _ = _ H (K E,A (r 2 )) _ 

det [diag([K e > A (Ti)] -1 ) diag(K E . A (T 2 ))] “ det[diag([K E . A (Ti)]-i)] • det(K E . A (Ti)) ’ 
det(R E ’ A (T 2 ,T!)) det(K E ’ A (T 2 )) 

det(diag(R E ’ A (T 2 ,T 1 ))) det(diag(R E - A (T 2 ,TO)) • det(K^ A (T!))’ 


where Eq. (l50l) is used. The subscript R denotes real, and subscript C denotes complex. Hr (R e ' A (T 2 , Ti)) is real 
because both its numerator and denominator are real. 


In order to quantitatively estimate the extent of resonance overlap, we use the same overlap matrix as in Ref. 
Q, but include only the |a) states, which are populated by the exciting laser spanning the energy range [El,Eh\: 

Ky= E l<«|5>H<a|K'>l- (53) 

a,E a E[EL,EH ] 

The Hadamard non-diagonality measure for the matrix of size Nq x Nq , composed of K , values, is introduced 
as 


H = det {fl a ) /det (diag (ft a )). 


(54) 


The numerator in Eq. (15411 is shown numerically to be always real and positive, and the denominator is equal to 
n,E K , and thus also real and positive. The same inequality as in Eq. (l49l) is valid for H (fi“). 


D. Implementation of the Shaped Laser as a Linear Combination of Gaussian Laser Pulses 


The eigenvector e A providing the desired optimized value A^ ,a after the pulse is over [Eq. (1461) ] is a finite 
discrete set of complex values of laser amplitudes £ p (u)A,g ), at different frequencies. These values can be reached in 
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multiple ways. The approach used for the IBr model m , is also used here: namely, to obtain the desired set of 
£ p (uiA,g ) values, A = 1,..., Na, it is sufficient to take the same number of linearly independent functions £ a (w), and 
expand the components of £ A in terms of £ a (w) at all tOA,g frequencies with the (as yet unknown) time-independent 
complex coefficients d a - 

£p[cO A,g ) — ^ ^ d a £ a {uJ A ,g ) , Q, A — 1,..., Na , t ^ Tover . (55) 

a 

or, as a matrix equation: 

£ — B d, PA,a £a(^ A,g) , d — {d 1 , . . . , dj\[ A ^ , t Tover■ (56) 

The set of £ a (w) functions is linearly independent, the B determinant is nonzero, and the unique nonzero vector d 
exists as a solution of Eq. found as 

d = [B] -1 £ A . (57) 

The basis functions £ Q (w) in frequency domain can be assumed to be infinite-time Fourier transforms of the corre¬ 
sponding basis functions £ a {t) in time domain (the latter are all vanishing when t > T over ). In turn, finite-time Fourier 
transforms of £ a (f) can be written as £ a (uj,t), and at finite times Eq. (l55l) takes the form: 

£p(^ L? J 4,p; t'j ^ ] d a £ a A,g, t) ‘ CL, A — 1,..., Na] (58) 

a 

i.e., 

£ A (f)=B(f)d, B A ,a(.t) = £ a (ujA,g,t), d = (di,.. . ,d NA ) T . (59) 

Using Eq. (15^1) . Pg 2 (f), Eq. (1751) . can be expressed in terms of the d vector: 

Ps 2 (t) = £ At (t)K e ’ A (t)£ A (t) = d t B t (t)K e ’ A (t)B(t)d. (60) 


Thus, the d vector in Eq. (1571) can be used for time propagation of Ps 2 ( t ) [Eq. (RJOl) ] at all times: before the 
laser is turned on, while the laser is on, and after the laser is off. Optimized populations always satisfy the condition 
Ps 2 (t = T 2 ) = \ R - A P S2 (t = T 1 ). 

To perform numerical computations, we select a set of Gaussian laser pulses £ a (t), centered at different frequen¬ 
cies uj a : 


£«(i) = e a /(2 V^cta) exp (- ( t /(2 a a )) 2 - ioJa^j ■ 

The finite-time Fourier transform of this Gaussian pulse, £ 0 (u;,f), [Eq. (1181) ]. can be expressed analytically 


(61) 


18-20] as: 


£a(w, t) = (e a /2) exp (w — w a ) 2 ^ |2 —exp (a a (uj - co a ) + it/(2a a )f fU(a a (w - u> a ) + zt/(2a 0 ))| , (62) 
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where W ( z ) is the complex error function 


[20, 21]. 


At times t > T over = 4V2 In 2 a a this becomes 


£a(w) = e a exp (-al (w - w a ) 2 ) . 


(63) 


Using Eq. (E3J, the control pulse £ P (t) in time domain is 

Na N a 

£ p (t ) = d a£a(t) = ^2 d a e a /(2V^a a ) exp (f/( 2a a )f - iu a tj , (64) 

a= 1 a— 1 

with infinite-time Fourier transform 

n a n a 

£ p {uj) = ^2 da,£a(u) = ^ d a C a exp (-al (w - Wa) 2 J ■ (65) 

a =1 a=1 

By construction, the £ p (uiA, g ) value should be constant inside the corresponding I a bin [Eq. (1 6 6 [) ]. The £ p (co) function 
[Eq. (1551) ] is smooth and does not satisfy this requirement exactly. Nevertheless, if Na is large enough each I a bin 
becomes relatively small, and the smooth function in Eq. © in each bin can be approximately treated as constant. 


IV. COMPUTATIONAL RESULTS 


Consider Sq -a S 2 excitation to coherently control S 2 tA Si interconversion dynamics of pyrazine excited using 
weak light in the perturbative regime. We use the pyrazine vibronic structure of Refs. [3] and ij], and partition the 
energy into 2000 bins, in the range 4.06-6.06 eV, where energy is referred to the ground vibrational Sq state. Here, 4.06 
eV is the Si energy at the So nuclear equilibrium configuration [12, 22]. The Q space consists of the 176 brightest (most 
optically accessible) |re) resonances, having the largest values of {n\ii\g). In this case the QP-partitioning approach 
gives 76775 coarse-grained vibronic states |cf), with energies ranging from 4.06 to 6.06 eV. Thus, there are 76775x176 
= 13512400 K = (a|re) values. These are used together with 176 (re|/z|g) values to compute the dynamics of interest. 


A. Uncontrolled Excitation and Decay Dynamics 


Figure |T] shows characteristic examples of Ps 2 (t) populations produced by a single Gaussian laser pulses of 
differing time durations, where the subscript u denotes “uncontrolled”. These examples are computed with the laser 
center frequency corresponding to 4.84 eV. It is notable that the uppermost population curve in Fig. [TJ produced by 
the pulse with a time duration ~1 fs ( a a = 0.1 fs) is, at times t > 0.5 fs, similar in shape to the zero-zero curve in 
Fig. 5, Ref. P|. This is the case because the ultrafast laser pulse behaves like e a S(t ) on the femtosecond timescale, 
and its finite-time Fourier transform is nearly constant, ~ e a . As a consequence, in this specific case, after the pulse 
is over, Ps 2 (t) in Eq. (1311) is the same up to a constant scaling factor as the zero-zero Ps 2 {t ) in Eq. (0, with 
C K > oc {i/h){n'\n\g)e a - 
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Figure [2] shows Ps 2 (t ) populations produced by Gaussian lasers having the same short time duration «10 fs 
( a a = 1.0 fs), but different center frequencies. In this case all populations behave similarly on a short time scale, 
differing by the overall magnitude due to the difference in (n'\n\g) values for different resonances | n'). 

Figure [3] shows Ps 2 (t) populations produced by Gaussian lasers with long time duration around 200 fs (a a = 20.0 
fs), using different frequencies. In contrast with Fig. [2j there are significant differences in S 2 Si IC dynamics, 
depending on the frequency used. Figure [3] shows that the laser with 4.84 eV photon energy produces a larger 
population, which also tends to decay slower, than in other cases, thus, marking the region of relative stability in 
pyrazine resonance structure. 


Ref. 


Both Figs. [2] and [3] qualitatively correlate well with the corresponding results for Sq —> S 2 O Si dynamics in 


23], obtained using a more general non-perturbative time-dependent dynamical approach 


241. 



FIG. 1: S 2 populations Ps 2 (t), denoted P u (t ) here, produced by Gaussian laser pulses of different time duration. Panel inset: 
The same data, shown on a shorter time scale. 
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FIG. 2: S 2 populations, Ps 2 {t), denoted P u (t) here, produced by short Gaussian laser pulses with the same a a = 1.0 fs, but 
different center frequencies. 


B. Control Involving Multiple Overlapping Resonances 

Consider first sample numerical results for H (fl“), the measure of the extent of resonance overlap [Eq. (1541) ] 
and the quantities associated with it. These quantities are H (K £,A (i)) [Eq. (1151) ]. which is the K £,A (t) non- 
diagonality measure, shown at Tf = 150 fs and T 2 = 250 fs; and two measures of the non-diagonality R £,a (T 2 ,Ti) , 
Hr (R £ ’ A (T 2 ,Ti)) [Eq. (l5ll) ]. and | He (R £,A (T 2 , Ti)) | [Eq. (1521) ]. In addition, we tabulate A^/^ and A^ A , denoting 
minimal and maximal eigenvalues of the eigenproblem in Eq. TO and which we term “control extents”. Values for 
128 I a bins (degrees of freedom of the laser), are listed in Table |T| Note first the enormous range of control possible 
for the ratio Ps 2 (T 2 ) / Ps 2 (Ti) as indicated by the A^j^ and A^ ,A . For example, for the first energy interval, this ratio 
can range from 3.05 x 10 -6 to 3.90 x 10 +5 , a range of over 1 x 10 +n . 

The measures in Table I are obtained using products of 128 matrix elements of the corresponding matrices. 
Since each of these values is small, we report the 1/128 power of these measures. From Table Q] one can see a 
well defined correlation between H (fi“) and the other quantities. Generally, when H (fl“) is small, so too are 
H (K £ ’ A (Ti)), H (K £,a (T 2 )) , Hr (R £,a (T 2 , Ti)) and | He (R £,A (T 2 , Ti)) | (meaning a larger extent of non-diagonality 
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FIG. 3: Bottom: S 2 populations, Ps 2 {t) > denoted P u {t) here, produced by long Gaussian laser pulses with the same a a = 20.0 
fs, but different center frequencies. 


TABLE I: The (1/128) power of H (fi“), H (K e ’ a (Ti)), H (K £ ’ a (T 2 )), Hr (R e ’ a (T 2 ,7\)), \H C (R £ ’ a (T 2 , Ti)) |, as well as 
A^. A and A r „’ A for different energy intervals [El, Eh]- Here, Tf = 150 fs, T 2 = 250 fs. 


[El, Eh], eV 

[4.46, 4.66] 

[4.66, 4.86] 

[4.86, 5.06] 

[5.06, 5.26] 

[5.26, 5.46] 

H(n a ) 

1.76x 10 _1 

2.77xl0 _1 

3.09X10” 1 

3.15xl0 _1 

2.80xl0 _1 

H (K e A (Ti)) 

1.08xl0 -2 

2.68xl0 -2 

9.05xl0~ 2 

1.36xl0 _1 

1.13xl0 _1 

R(K eA (T 2 )) 

1.23x 10 -2 

2.50xl0 -2 

9.45xl0" 2 

1.29X10" 1 

1.06X10" 1 

Hr( R e ’ a (T 2 ,Ti)) 

1.41 x 10~ 4 

8.00X10" 4 

1.20x 10~ 2 

2.38xl0“ 2 

1.51xl0“ 2 

1 He (R e ' a (T 2 ,Ti))| 

1.36xl0 -4 

9.05X10" 4 

1.75xl0" 2 

3.41xl0“ 2 

1.79xl0“ 2 

\R,A 

A min 

3.05 x 10 -6 

3.36 xlO -5 

5.54x 10” 4 

1.29xl0“ 3 

7.30xl0“ 4 

\ R , A 

^max 

3.90x10+® 4.32xl0+ 4 

1.89xl0+ 3 

6.67X10+ 2 

1.92xl0+ 3 


in the corresponding matrices). In particular, correlation is good with A^ ,A — A r ^' irl ; when it is large, a greater extent 
of coherent control is possible, in agreement with the non-diagonality measures. 

Numerically implementing controlled Ps 2 (t) dynamics proceeded as follows. First, the eigenvalue problem in 
Eq. (1461) is numerically solved for the particular number of bins Na in the desired energy range [El, Eh], providing the 
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set of eigenvalues \ R A and corresponding eigenvectors e A , which give the X R ’ A as Ps 2 (T 2 )/Ps 2 (Ti) ratios during the 
P Sa (t) time propagation (P 2 > T\ > T over ). Then, a set of linearly independent Na Gaussian lasers [Eqs. (lUTI) (1651) ]. 
is introduced (all with the same a a ), contiguously and uniformly covering the desired energy range [El, Eh]- The 
eigenvectors obtained e A are then expanded in terms of this Gaussian basis with the d coefficients given by Eq. <53- 
The dynamics are then propagated from t < —T over to t > T 2 using the corresponding d coefficients for each e A 
eigenvector; finite-time Fourier transforms of the pulses in Eq. (1581) are produced using Eq. m, and the pulse time 
profiles are given in Eq. (PI) . The perturbative nature of the dynamics makes it possible to scale Pg 2 (t) uniformly by 
multiplying the e A eigenvector by a scalar constant. We utilized this scaling option to allow presentation of both the 
maximization and minimization results to be shown on the same figure (upper panel, Fig. |5j) below. Specifically, the 
maximization curve is multiplied throughout by 8.3 x 10” 5 

An experimental suggestion of R. J. Gordon (University of Illinois, Chicago) prompted our using a controllable 
laser in the wavelength range 250-265 nm, with time duration ~150-200 fs, to study the pyrazine So —»• iS 2 -O- Si 
excitation and IC dynamics. Using this as a guide, we computed control and dynamics in the corresponding energy 
range (El = 4.68 eV, Eh = 4.96 eV), using Tf = 150 fs, T 2 = 250 fs, Na = 128, and all a a = 21.0 fs. The resulting S 2 
populations, together with resulting control fields in time domain, are shown in Fig. 3] where the subscript c denotes 
“controlled”. Corresponding control fields in the frequency domain are shown in Figs. [5] and [6] 

The behavior of the controlled Ps 2 (t) (Fig. [1}, differs in magnitude in the regions when the pulse is acting, and 
after the pulse is over. To understand this difference, note that to obtain the controlled fields in Figs. [5] and [5] using 
a set of Gaussians requires that some components of d vector be large. After the pulse is over, these components 
are “balanced” by one another in the infinite-time Fourier transform, to give the small desired population value 
Pq at t = T\ or t = T 2 and to yield the required controlled dynamics. However, while the pulse is acting, these 
components are “unbalanced” giving large transient s p (u>,t) values. For similar reasons the controlled pulses, being a 
linear combinations of single Gaussians, are effectively longer than the single Gaussian pulse (see Fig. [I] lower panel). 

To examine the complex structure of the control pulses at Figs. [5]and[6l we apply several approaches to simplify 
the field while monitoring the control achieved. First, we attempted a local averaging of the controlled field, where 
the total field in Na bins is arithmetically averaged (amplitude and phase separately) using a smaller number Ns of 
larger bins (Na being an integer multiple of Ns, for example, for Na = 64, Ns = 32, 16, 8, 4, 2). By doing so, the 
resulting averaged field, however, showed virtually no control. Second, this averaged step-like field was expanded with 
Ns Gaussians and the resulting smoothed field used for the propagation. Again, this case led to nearly complete loss 
of control. 

An alternative simplifying approach was, however, successful. Specifically, we retained only the Nr largest field 
amplitudes out of the total Na (with all the smaller ampitudes set to zero), keeping the phase profile intact, and 
monitoring the changes in control ratios. Sample results for Na = 64 are shown in Fig. 7. A total Na of 64 is used 
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FIG. 4: Upper panel: Two controlled S 2 populations, Ps 2 (t), denoted P c (t), which either minimize and maximize \ R ’ A , i.e., the 
S 2 population ratio at times T 2 = 250 fs and T\ = 150 fs. The p™ ax curve has been multiplied by 8.3 x 10 -5 in order to fit on 
this figure. Lower panel: Time envelopes of two corresponding controlled laser pulses, |e p (t)|, together with the time envelope 
of the single (uncontrolled) Gaussian laser pulse, |e a (t)|. 

here (results with Na = 128 are qualitatively the same). It is clear from Fig. [7| that this approach, retaining only 
the largest amplitudes, works better than the previous two since it tends to partially maintain important dynamical 
information. Generally, A^)^ is more robust with respect to this amplitude truncation than is A^j^. Additionally, 
we found that the extent of control achieved using only Nr amplitudes out of Na, is similar in magnitude to control 
extents without truncation, but using this Nr as the original Na- That is, the same number of degrees of freedom in 
both cases provides similar extents of control. 

Theoretically, maximum and minimum control limits via this approach can be reached using all coarse-grained 
|a) states accessible to the laser, i.e., those belonging to the interval of interest [El, Eh]- For the case presented in 
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Amplitude profile of minimizing pulse 
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FIG. 5: Amplitude and phase of eigenvector, which minimizes the Ps 2 {T 2 )/Ps 2 (Ti) ratio. = 8.28x10 5 . 


Figs. HIE] and El the number of |a) states, using our pyrazine description, is 11885. However, as mentioned in the 
Appendix, the optimization problem for |cf) states in Eq. (144|) is numerically stable only up to dimensionality 150-180, 
and the control range continues to increase when the dimensionality increases from 128 to 180, reaching 

~ 10 5 . We anticipate a theoretical control range limit to be ~ 10 9 -10 10 , which, however, is not achieved due to the 
numerical limitations discussed in Appendix (see below). 


V. SUMMARY AND CONCLUSIONS 


Coherent control of internal conversion (IC) between the first and second singlet excited electronic states of 
pyrazine (Si and S 2 ) is examined, using two different control objectives. The control is performed by means of 
shaping the laser, which excites the system from the ground electronic state So to the second excited electronic state 
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Amplitude profile of maximizing pulse 
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FIG. 6: Amplitude and phase of eigenvector, which maximizes the Ps 2 (T2) / Ps 2 (Ti) ratio. A^’a^ = 8.38xl0 +3 . 


*S *2 • Resonance energy broadening and resonance overlap are shown to be responsible for phase control efficiency, and a 
correlation between resonance overlap and controllability is established. A huge range of control was obtained for the 
relative population of S 2 at long times as compared to times just after the pulse is over. Different ways to simplify the 
controlled fields are described, and the behavior of the control as a consequence of these simplifications is investigated. 
Specifically, we have found that retaining the largest field amplitudes is the best approach to field simplification. 
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FIG. 7: Upper panel: Dependence of A^j^ on the number of retained amplitudes Nr. Lower panel: The same, but for 
Total number of amplitudes Na = 64. 
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Appendix: Energy Partitioning for Control Associated with Overlapping Resonances 

As noted in Sect. IIII Bl numerical instabilities necessitated that we introduce a further partitioning of the 
energy axis. Specifically, we partitioned the energy axis into a limited number of Na bins in Eq. m- 


A ckG/a 


' 'y v £p{uA,gi t) 


E 

,a£l A 


= ^2 S p(^A,g,t)M e K ’A(t), 


where I a is a bin number A, which has the center energy Ea , u >A,g = ( Ea — E g )/h , and 

K^it) = E M^(t) 

aelA 


( 66 ) 


(67) 


is the collective material system matrix element, corresponding to bin I a- Here and below, the bold superscript A 
denotes that all the corresponding quantities are written for |A) states, which are defined below. 

Using Eq. (|28)l . M^'^(t) can be written as: 


= E («|a)r„(t)^E< a i ,s '>< K, i^> = ( K \ 

o;G/a «' 


E r a (t)\a){a\ 

.ckG/a 




( 68 ) 


The middle expression in square brackets, unlike Eq. (051) . is not the single |a) state propagator, but the localized 
coarse-grained propagator, with the sum only over |a) states belonging to the bin I a ■ 

One can introduce the “binned” states |A), such that the corresponding projector onto the state |A) is 

\a)(a\ = (i/n I a ) E I s )(“I> hence V^\ a )(M\/n~u = E 

q;G/a ckE/a 

where Nj a is the number of |a) states that are inside bin I a- With the notation |A) = ^ Ni a \ A) we have 


\A){A\ = E !«)(«!• 

oeU 

Using Eq. (l69l) . the propagator in Eq. (f68li can be approximately rewritten as: 

E r a (t)\a)(a\ « |A)<A| E *<*$) = r^(t)\A)(A\, 


(69) 


(70) 


ae Ia «G/a 

where r^it) = (1/Ni a ) J2 a ei A T a(t)- The accuracy of the approximation made in Eq. (l70l) rapidly increases with 
decreasing bin size. Using Eq. (1701) . M f E(t) in Eq. (l68l) can be rewritten as: 


~ {k\a)t% (t)^E( A i K ')<«Vl5>- 


(71) 


Defining 


9a = 1 E<4«'X* , M3> = t(A\9\9), R%k = (A\k) 


(72) 
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gives 

M E-A (f) = R At Z A (i)/A K E ’ A (f) = M E ’ At (t)M E ’ A (i) = /£ At z At (i)R A R. At Z A (t)/A (73) 

where r A (f) is a square diagonal matrix composed of r A (t) values, and /i A is a square diagonal matrix composed of 
(i/Tb)(A\n\g) values. This gives the Ps 2 {t) population in terms of binned values as 

Ps 2 (t) = e At (% At z At (f)R A R At z A (t) t£ A e A {t) (74) 

= , y \ I £p(bJA',gi t ) I K A , A i{t) + £p(a )A',g,t)e p (uJA",gP)pA' (75) 

A' A'^A" 

where £ A (f) is a vector composed of £ p (u>A,g,t) components. 

Since /r A and r A (t) are diagonal, the only possible source of nondiagonality in Eq. (1731) for K E,A (f) and Eq. 
d is Q a = R A R A t, composed of Q A , A „ = (A \Q\A ) values. Thus, the possibility of phase control by means of 
phases (f>A{t) of complex £ p (u)A,g,t) = \^p(^>A,git)\ exp (icf>A{t)) depends solely on its properties. As in the previous case 
of 1 7 ) and |a), all the Ps 2 {t) phase control considerations remain the same, except that |y) or |a) states are replaced by 
| A) states. Namely, phase control is provided by resonance energy broadening and resonance overlap. The resonance 
overlap effect, providing the non-block-diagional structure of Q A and K e,A (f) as a consequence, enhances the effect 
of resonance broadening. 

Using Eq. dm the eigenproblem in Eq. (l44| is reformulated as 

R e ' A (T 2 ,Ti)£ A = \ R ’ A e a , 

R e ' A (T 2 ,Ti) = [K e,a (T 1 )] _ 1 K e,a (T 2 ). 


(76) 

(77) 


Its dimensionality reduced from N a to Na , allowing an accurate numerical solution for Na values up to 150-180. 


